Fluctuation-Induced Casimir Forces in Granular Fluids. 
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We have numerically investigated the behavior of driven non-cohesive granular media and found 
that two fixed large intruder particles, immersed in a sea of small particles, experience, in addition 
to a short range depletion force, a long range repulsive force. The observed long range interaction 
is fluctuation-induced and we propose a mechanism similar to the Casimir effect that generates 
it: the hydrodynamic fluctuations are geometrically conflned between the intruders, producing an 
unbalanced renormalized pressure. An estimation based on computing the possible Fourier modes 
explains the repulsive force and is in qualitative agreement with the simulations. 
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Granular materials have been extensively investigated 
because of their complex dynamics [l| . Examples of this 
include: pattern formation, Faraday waves, avalanches, 
convection phenomena, segregation and many more. One 
of the most active fields is the behavior of granular mix- 
tures, e.g. Brazil nut effect J3| and its multiple varia- 
tions [3|. In these experiments, the granular material is 
composed of a mixture of two types of particles differ- 
ing in mass or size. When the system is agitated, parti- 
cles of different types may group together (demixing) or 
stay mixed [^ . Phase diagrams for the mixing/demixing 
transition have been constructed for different material 
properties or experimental conditions. However, a fun- 
damental question remains unanswered: Is the mixing or 
demixing caused by an effective long range force between 
the particles? If so, which is its origin? Let us note that 
a long range force is difficult to justify a priori, as grains 
interact only through a short ranged hard-core potential. 

Recently, three experiments on driven granular mix- 
tures |5|, la, I3 1 performed under very different experimen- 
tal conditions, give a hint on how to answer the previous 
question. These experiments have shown that thermody- 
namic properties (like pressure |a| , density |^ and veloc- 
ity fluctuations [TJ) are different in the regions between 
the larger particles, versus the remaining, external, re- 
gions. Consequently, the big particles modify some phys- 
ical properties in the confined area between them. The 
most likely reason is that the larger particles limit the 
allowed wavevectors of the hydrodynamic fluctuations of 
the small particles that surround the larger ones. 

Forces arising from the conflnement of a fluctuation 
spectrum have attracted attention since the seminal work 
of Casimir in 1948, who predicted the existence of an 
attractive force between two metal plates, separated by 
a vacuum, due to constraints on the quantum electro- 
magnetic field in the gap imposed by the conducting 
plates y|. In fact, the concept of Casimir force is 



more general and is common to systems characterized 
by long-range fluctuations subject to a geometrical con- 
straint which limits the long- wavelength portion of their 
spectrum. Soft condensed matter provides examples of 
Casimir forces, such us those arising in confined criti- 
cal fluids, in liquid crystals and superconducting films, 
where long-range correlations are the consequence of a 
broken continuous symmetry |2j . Vibrated granular mat- 
ter ll(l | and granular avalanches [llj provide other exam- 
ples of physical system where correlations can become 
long ranged, in spite of having short-range forces. 

In this Letter we show that there is an effective long 
range force between the large and heavy particles in a 
granular mixture. Two ingredients are required: (i) long 
range correlations and (ii) the confinement of the fluc- 
tuation spectrum induced by the large particles in the 
density, velocity and temperature fields. 

We consider the driven granular model in [lO|, Il2l |. 
Grains are hard particles of diameter d and mass m and 
their collisions are characterized by a constant normal 
restitution coefficient a. To achieve a stationary state, 
energy is supplied into the system by random forces act- 
ing on all particles. The random forces F, are modeled as 
a white noise of intensity F: (Fi(i)Ffe(i')) = mT5ik5{t — 
t'). This system reaches a homogeneous stationary state 
characterized by long range correlations, leading to the 
renormalization of the energy density and collision fre- 
quency due to the fluctuations at low wavevectors [l(j|. 

The system is composed of N grains put in a square 
box with periodic boundary conditions. Besides the small 
grains, two inelastic impenetrable and immobile large 
hard disks (the intruders) of diameter D are placed, sep- 
arated at a distance R. The coefficient of restitution 
a is the same for all types of collisions. This granu- 
lar mixture is studied using molecular dynamics simula- 
tions. The grain-grain and grain- intruder collisions are 
treated as usual, using an event-driven code. To take 



into account the random forces, a new type of event is 
introduced: collisions between particles and a "thermal 
bath" . The event results in a momentum p being in- 
stantaneously transferred to the particle, where p is ran- 
domly chosen by sampling a Gaussian distribution with 
zero mean and given variance -P^ath ^^ each direction. In- 
teractions with the thermal bath are scheduled for each 
grain. When an interaction takes place for a particle, 
its momentum gets updated and a new future event is 
scheduled for the same particle, after a random interval 
of time tnoxt, according to an exponential distribution 
-P(incxt) ~ exp(-inoxt/Tbath)- In the hmit Tbath -^ and 
^bath — > 0, this injection method converges to the white 
noise force with T = -Pbath/"^''bath- In practice, the time- 
scale for interaction with the thermal bath, Tbath, is taken 
smaller than the frec-flight-time. 

Hereafter, we choose as basic units d, m, and F. These 
units define the time unit as to = {mcP/TY^^ and en- 
ergy unit as cq — {mcPT^y/'^. We take the diameter 
of the intruders D = 8d and the coefficient of restitu- 
tion a = 0.8. We simulate systems of size L = 60d and 
L — 80d, with the number density of the granular fluid 
n — 0.366 d~^. Given the density, restitution coefficient 
and noise intensity, the stationary temperature can be 
computed using mean field models giving Tq = 1.84 cq, 
and the collision frequency is vq — 3.03 1,^ . Hydrody- 
namic fluctuations determine a stationary temperature 
higher than Tq that depends on the system size [ifl • For 
L ^ 60d, T = 2.43 eo and for L = 80d, T = 2.46 cq. 
For every configuration, simulations were run for about 
5 X 10^ collisions per particle. We investigated the ef- 
fective interaction between the intruders at a distance 
R. For this purpose, we measured the component of the 
total momentum transferred from the gas to intruders 1 
and 2 along the line, parallel to the x-axis, joining their 
centers, Pix {i = 1,2) as an average over a time inter- 
val T — 33.4 to (corresponding to about 100 colhsions per 
particle). This procedure gives the "instantaneous" value 
of the fluctuating force as F12 = {P2X — Pix)/'2t, whose 
time-average finally leads to the net effective force, F. In 
the elastic case, a = l,r — 0, F vanishes as expected, 
whereas for a < l,r 7^ 0, -F is definitively different from 
zero, showing an effective force between intruders. The 
average y-component of the force is compatible with zero 
to numerical accuracy. Varying the noise strength, we 
checked that the force is proportional to the granular 
temperature, as can be deduced by dimensional analysis. 

Figure ^ shows the relative force as a function of dis- 
tance R for the two system sizes. The inset shows the 
short range part of the relative force, that extends for 
some small particle diameters, d, alternating between at- 
traction and repulsion. At larger distances, a repulsive 
force is observed with an interaction range much larger 
than d and comparable with D or the box width (actu- 
ally, due to the periodic boundary conditions, the force 
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FIG. 1: Relative effective dimensionless force FD/Tq be- 
tween the large particles as a function of the distance R be- 
tween their centers. We plot simulations results for L = 60d 
(open triangles) and for L = 80d (filled circles), together with 
theoretical predictions (solid lines). Inset: short range region 
showing the depletion forces. Note the difference in the ver- 
tical scale, compared with the long-range part. Error bars in 
the inset are smaller than the symbols. 



must vanish at i? = L/2, as seen from the simulations). 

The oscillatory short range part is similar to the de- 
pletion forces appearing in fluids that develop local lay- 
ering structures [ij, [ij], which are usually explained 
by entropic arguments based on equilibrium statistical 
mechanics. Recent works 6, 15] suggest that depletion 
forces might be at work in granular mixtures, so the con- 
cept of entropic forces is applicable here. However, the 
long range part cannot be due to a depletion mechanism 
as it extends beyond the typical range of the depletion 
forces, and because it is repulsive in all its range. 

To elucidate the nature of this force, we analyzed the 
probability distribution of the force over the intruders, 
plotted in Fig. [3 There it is seen that fluctuations are 
about 20 times larger than the average force, and there- 
fore very long simulations are required. As the force 
is proportional to the granular temperature, these large 
fluctuations are not an artifact of a high temperature, 
but an intrinsic property. Large fluctuations is a generic 
signature of fluctuation- induced forces (see, e.g., |2l|la|)- 

This fact together with the known property of large 
fluctuations in driven granular media, suggest that the 
long range repulsion is a fluctuation-induced force as in 
the Casimir effect [sl |9j. The confinement of the hy- 
drodynamic fluctuations between intruders restricts the 
allowed fluctuation modes to those with wavelengths 
smaller than the gap size, whereas the spectrum of fluc- 
tuations in the outer region allows smaller wavevectors 
and forms a quasi-continuum, as illustrated in Fig. |31 
Consequently, the hydrodynamically-generated "radia- 
tion pressure" between the intruders is different from the 
one outside this inner region. 

To describe the Casimir force originated by the flue- 
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FIG. 2: Probability distribution of the fluctuating force F\2 
aX R = 20d, for L = 80d. The net effective force is obtained 
as the average of this histogram. The average value is F = 
0.223To/D and the standard deviation is ap = 4.83ro/D, 
that is about 20 larger than the average. 



tuating hydrodynamic fields on the intruders we use an 
approach similar to Ref. [lOJ , where the pressure is renor- 
malized in each point due to fluctuations that are com- 
puted using fluctuating hydrodynamics. Instantaneously, 
the pressure tensor at position r is given by 

P*{y) ^p{n{r),T{r))I + mn{r) u(r) u(r), (1) 

where n and T, and u are the fluctuating density, tem- 
perature, and velocity fields and / the identity ten- 
sor, respectively. Moreover, p{n,T) = TH(n) is the 
usual thermodynamic pressure for hard disks [13 with 
H{n) = n(l + 02/8)/(l - c/))^ , and cj) = 7m(f/A. As the 
intruders are immobile, u vanishes at their surface, so 
the contribution of the convective term in p* vanishes 
and it becomes a scalar. This would not be the case if 
the intruders were allowed to move. 

Linearizing the hydrodynamic fields (n, T, u) around 
the stationary values, (ng, Tq, 0), we can expand the pres- 
sure up to second order in the fluctuations Sn and ST. 
Its statistical average over the random noise is 



{p*) - po + HiiSTSn) + ToH2{Sn^), 



(2) 



where 



Po 
1^21 



p{no,To), Hi 



dH/dn\ng, and 



H2 — ^d H/dn Ing. The density-density and density- 
temperature fluctuations appearing in Q are difficult to 
compute because it is needed to solve the fiuctuating hy- 
drodynamic equations with the full boundary conditions 
imposed by the intruders. Alternatively, one could solve 
the corresponding equation for the correlation functions 
including the boundary conditions. A simpler estimation 
of (p*) can be obtained by employing the Fourier trans- 
forms of the fiuctuating fields M(r) — V~^ J2k e~''^'' ^SAu 
and the structure factors SabO^) = V^^{SAis^ (©-k) |lO|. 
Expression ^ transforms into: 



(P*> =Po + V-^ Y. [^^ ^"^('^) + ^0 ^2 ^™(k)] . (3) 




FIG. 3: Sketch of the hydrodynamic fluctuations leading to 
the Casimir-like force. In the gap between the intruders only 
fluctuations of wavelength smaller than R — D are allowed, 
while in the outer space the wavelengths can be larger. 



The structure factors for the uniformly driven system 
have been described in detail in [lOj. The relevant con- 
tribution comes from the region at small k, where they 
show a power law dependence SabO^) = S^g k~^. The 
prefactors, 5^^, depend on density, noise intensity, and 
restitution coefficient a. These asymptotic expressions, 
when inserted into ^, yield a sum C^j^. fc^^, where the 
coefficient C = Hi S'JJy -f Tq H2 5'^„ turns out to be neg- 
ative for n < 0.73 d" 2 and positive for larger densities. 
Therefore, for n < 0.73 d"^ (the case of the simulations), 
fluctuations produce a decrease of the local pressure. In 
the gap between the intruders (region I according to Fig. 
|3J) the number of possible k modes of low k is smaller 
than outside (region II). The effect is that the pressure is 
lower outside than inside, leading to an effective repulsive 
force between the intruders. Furthermore, C is propor- 
tional to the temperature, so is the Casimir force. The 
sign of the force is reversed for densities n > 0.73 d"^, 
which are too close to the freezing transition or random 
close packing [l8| to be observed in the simulations. Ex- 
perimentally, analogous crossover is found by increasing 
the driving intensity I3. 

Due to the long range correlations, fluctuations in 
regions I and II of Fig. O are correlated. However, 
in order to numerically estimate the /c-sums, we treat 
these regions as being independent. This approximation 
will overestimate the pressure difference and hence the 
Casimir force. However it provides a rough estimation of 
its numerical value. In detail, we perform the fc-sum only 
over the fc-vectors allowed by the geometrical constraints. 
In a rectangular box of size axb, the x-component of the 
k vectors is 2Trnx/a and the y component is 2Tmy/b. It is 
at this point where the difference between regions I and 
II appears: a = R — D m region I and a = L — R — D 
in region II, while 6 = _D in both regions. Moreover, the 
vector k = (0, 0) must be excluded from the sum, and we 
introduce an ultraviolet cutoff, kc, beyond which hydro- 
dynamics is not valid. We take the cutoff kc — 2TT/do, 
where do — max((i, Zq), and Iq is the mean free path of 
the small particles. 

In the limit of small fc, with structure factors going 



as A;~^, the pressure (O can be analyzed asymptotically 
in the cases a ^ b (for particles at large distances) and 
a <C 6 (for particles at short distances) 



{P*) = 



Po + C a/b ; a 3> fe > do 
Po + Cb/a ; 6 3> a > do. 



(4) 



Note that these asymptotic expressions for the renormal- 
ized pressure do not depend on the cutoff distance, as 
long as a and b are much larger than it. Finally, the effec- 
tive force on the particle 2 is the difference of the forces at 
the left and the right of the particle F2 = D[{p^) - (pii)]- 
A negative value of C gives rise to a long range linear 
repulsive force and a short range attractive force, at dis- 
tances smaller than D. The opposite is obtained when C 
is positive. Note that in this estimate the force depends 
on the system size. This fact is related to the structure 
of the fluctuations, that become larger for small wavevec- 
tors. Therefore, increasing the system size, while keeping 
R fixed, the fluctuations in region II become larger, de- 
creasing even more the pressure in this region. However, 
as mentioned, at long distances this approximation is not 
completely valid. 

The Casimir force can be computed numerically us- 
ing the full expression of the structure factors [l(J, not 
only the k~^ part. The computed force is shown in Fig.^ 
Note that the linear force dependence is preserved but the 
attractive part is lost, mainly because the simulation box 
is finite and corrections to the fc~^ order are observed. 
The computation using the k~^ part shows a small dif- 
ference of about 10% compared to the full calculation. 
The comparison with the simulations indicates that this 
prediction overestimates, as expected, the Casimir force, 
especially in the L — 80d case, but gives the same order 
of magnitude and correct sign of the force. Moreover the 
predicted force for L = 80d is larger that for L ~ 60d, 
in agreement with the simulations. However, our the- 
ory does not predict the saturation of the force observed 
for R > 20d. Possible sources for this discrepancy are: 
(a) hydrodynamic correlations between regions I and II 
reducing the pressure difference, (b) geometrical factors 
that arise from considering rectangular regions instead of 
those bounded by circles, and (c) fixed intruders break 
the Galilean invariance and may modify the structure 
factors at very short wavelengths. 

To summarize we have found that two intruders, im- 
mersed in a sea of smaller granular particles driven by a 
white noise force, experience a long range mutual repul- 
sion. This repulsion has a dynamical origin and cannot 
be explained by standard depletion forces. We have pro- 
posed a mechanism based on the confinement of hydro- 
dynamic fluctuations when the intruders are near. The 
present effect is an example of repulsion determined by 
fluctuation-induced forces instead of the standard attrac- 
tion; a phenomenon which has been predicted to occur 
also in polymers [l9j. 



We claim that the force we observe is the granular ana- 
log of the Casimir effect. We propose a novel method, 
valid for non-equilibrium fluids, to compute Casimir 
forces starting from the structure factors. The two key in- 
gredients which render this effect manifest in the context 
of granular gases are: (a) the occurrence of large low-fc 
fluctuations, originating in the coupling of nonconserving 
noise with conserving hydrodynamic fields (conservation 
of particle number and momentum in collisions) [10|| , (b) 
the confinement of these fiuctuations in a gap (the space 
between the intruders), which is the common feature of 
all instances of Casimir forces. The enhancement of the 
low-fc components plays the role of criticality in equilib- 
rium molecular systems, where forces are induced by the 
thermal fluctuations of a correlated fluid in a confining 
geometry. Finally, these long-range forces might be re- 
sponsible for segregation effects in vibrofluidized granular 
mixtures of particles having different material properties. 
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